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考虑 客户 价值 的 卡车 与 无 人 机 联合 配送 时 变 路 径 优化 方法 研究 
温 廷 新 ， 吕 艳 华 ; 


(辽宁 工程 技术 大 学 工商 管理 学 院 , 辽宁 戎 芦 岛 125000) 


摘 要 : 针对 车 辆 拥有 量 的 增加 和 拥堵 情况 的 日 益 严 重 导 致 的 物流 业 配送 时 效 不 高 、 客 户 价值 低 等 问题 ， 综 合 考虑 
客户 价值 和 成 本 等 因素 ， 提 出 了 一 种 卡车 与 无 人 机 联合 配送 时 变 路 径 的 优化 方法 。 考 虑 到 配送 过 程 中 不 同时 段 的 拥 
堵 情 况 , 采用 速度 分 布 函 数 刻 画 车 辆 的 行驶 速度 , 同时 考虑 客户 的 时 间 窗 、 车 辆 的 载重 和 无 人 机 的 载重 等 约束 条 件 ， 
建立 了 成 本 最 小 的 数学 模型 。 根 据 模 型 的 特点 ， 引 入 k-means 对 客户 的 位 置 进行 聚 类 ， 设 计 混 合 的 粒子 群 算法 对 模 
型 进行 求解 。 最 后 通过 Solomom 数据 进行 模拟 仿真 实验 ， 对 模型 和 算法 的 有 效 性 进行 验证 。 实 验 结果 表明 ， 与 未 考 
虑 客户 价值 静态 路 网 模型 相 比 ， 该 模型 在 降低 9.32% 成 本 的 情况 下 ， 同 时 提高 了 16.83% 的 客户 价值 和 21.28% 的 客 
户 满意 度 ， 所 提 算 法 在 降低 配送 成 本 和 提高 企业 经 济 效益 方面 具有 一 定 的 有 效 性 。 

关键 词 : 车 载 无 人 机 ; 客户 价值 ; 时 变 路 网 ; k-means 聚 类 ; 联合 配送 ; 混合 粒子 群 算法 
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Research on time-varying route optimization method for truck and UAV joint delivery considering customer value 


i Wen Tingxin, Lyu Yanhuai 
(School of Business Administration, Liaoning Technical University, Huludao Liaoning 125000, China) 


Abstract: Aiming at the problems of low timelessness and low customer value caused by the increase of vehicle ownership 
and congestion, the paper proposes a time-varying route optimization method of truck and UAYV joint delivery after considering 
the factors of customer value and cost. Considering the congestion in different time periods in the delivery process, it uses the 
speed distribution function to describe the vehicle speed and the constraints such as customer time window and it takes vehicle 
load and UAYV load into account to establish a mathematical model with minimum cost. According to the characteristics of the 
model, it introduces k-means to cluster customers' locations and designes a hybrid particle swarm optimization algorithm to 
solve the model. Finally, the paper uses Solomom data to conduct simulation experiments to verify the validity of the model 
and algorithm. The experimental results show that compared with static network model without considering customer value, 
the model can reduce 9.32% cost, improve 16.83% customer value and 21.28% customer satisfaction, and the algorithm 


proposed in this paper has certain effectiveness in reducing distribution cost and improving enterprise economic benefits. 
Key words: vehicle UAV; customer value; time-varying road network; k- means clustering; joint distribution; hybrid particle 
swarm optimization 


下 ， 将 多 车 场 作为 配送 的 场景 


六 | 


系 。 


到 局 为 进一步 求解 数学 模型 ， 多 数学 者 采用 改进 的 启发 式 智 
随 着 电子 商务 的 发 展 ， 如 何 提高 商品 配送 效率 、 降 低 物 ”能 优化 算法 进行 求解 。Salama 等 四 为 优化 卡车 与 无 人 机 的 配 

流 成 本 ， 成 为 整个 物流 行业 面临 的 关键 问题 。 在 此 背景 下 ， 送 路 线 ， 提 出 了 一 种 基于 无 监督 机 器 学 习 的 启发 式 算法 来 加 

技术 的 发 展 和 需求 的 更 迭 ， 将 高 效率 、 低 成 本 的 无 人 机 用 于 ”” 快 求解 速度 。 王 新 等 入 为 提高 客户 的 满意 度 ， 综 合 考虑 无 人 

物流 配送 成 为 物流 行业 内 研究 的 热点 。 然 而 ， 无 人 机 存在 负 ”机 站 点 和 客户 时 间 窗 要 求 ， 建 立 以 总 成 本 最 小 化 为 目标 的 问 


载 小 、 飞 行距 离 短 等 局 限 ， 导 致 无 人 机 无 法 独立 完成 大 规模 。 题 模型 ， 并 设计 自 适应 大 规模 邻 域 搜索 算法 进行 求解 。 邓 水 
的 物流 配送 任务 ， 因 此 卡车 与 无 人 机 联合 配送 模式 成 为 物流 。 获 外 等 在 自然 灾害 情境 下 建立 卡车 与 无 人 机 联合 配送 冷 链 物 
行业 降 本 增 效 的 必然 选择 。 流 优化 模型 ， 采 用 进化 逆转 操作 ， 并 设计 改进 的 遗传 算法 。 

近年 来 ， 国 内 外 关于 卡车 与 无 人 机 联合 配送 路 径 问题 的 上 唐 立 等 (WI 利 用 泰 森 多 边 物理 概念 ， 改 进 蚁 群 算法 ， 对 考虑 路 
研究 主要 针对 模型 的 建立 和 求解 算法 的 设计 两 方面 。 从 模型 ” 径 安 全 的 山区 无 人 机 路 径 模型 进行 求解 。 李 妍 峰 等 1 改进 变 


的 建立 方面 ，Murray 等 由 首次 将 卡车 和 无 人 机 结合 ， 建 立 了 邻 域 搜索 算法 求解 需求 可 拆 分 的 路 径 问 题 。 曹 英 英 等 9 利用 
卡车 与 无 人 机 联合 的 旅行 商 问题 的 混合 整数 规划 数学 模型 。 遗传 模拟 退火 两 阶段 算法 求解 集群 下 的 卡车 与 无 人 机 联合 配 
在 此 基础 上 , Semiz 等 所 将 客户 时 间 窗 作为 约束 条 件 , 建立 了 送 问 题 。 能 兴隆 等 0 以 配送 时 间 最 短 为 目标 ， 分 为 两 步 ， 提 
带 时 间 窗 的 卡车 与 无 人 机 的 路 径 模 型 。 朱 晓 宁 等 中 根据 车 辆 出 新 型 优化 迭代 算法 进行 路 线 的 规划 。Han 等 风采 用 一 种 改 


限行 和 空域 禁 飞 的 情况 ， 将 区 域 限制 因素 嵌入 到 模型 的 构建 ” 进 的 人 工蜂 群 算法 来 最 小 化 系统 的 整体 运行 成 本 。 Phan 等 03 
当中 。 茧 勇 等 和 由 在 疫情 背景 下 ， 以 配送 商品 时 间 最 短 为 优化 。 将 问题 扩展 到 一 辆 卡车 与 多 架 无 人 机 结合 的 情形 ， 采 用 贪 梦 
标 ， 设 计 混合 邻 域 搜索 算法 求解 无 人 机 为 多 个 客户 无 接触 ” 随机 自 适 应 搜索 程序 (GRASP) 和 自 适应 大 邻 域 搜 索 (ALNS) 
送 的 路 径 问题 。 张 得 志 等 回 从 低 碳 与 随机 需求 两 方面 出 发 ， ”启发 式 算法 来 解决 这 一 问题 。 以 上 研究 从 成 本 最 低 和 距离 
立 了 多 目标 模型 。 杜 茂 康 等 四 在 实际 物流 配送 要 求 的 前 提 ” ” 短 目 标 出 发 ， 但 均 未 考虑 客户 价值 因素 ， 同 时 ， 上 述 研究 
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录用 定稿 


是 在 静态 路 网 环境 下 对 路 径 优 化 的 探索 ， 一 定 程度 上 忽略 了 
时 变 路 网 对 物流 配送 优化 的 影响 ， 建 立 的 模型 鲁 棒 性 较 差 ， 
会 增加 物流 运输 企业 的 配送 成 本 。 

基于 此 ， 本 文 在 满足 车 辆 载重 、 客 户 时 间 窗 、 无 人 机 飞行 
距离 和 无 人 机 载重 的 约束 条 件 下 ， 将 客户 价值 作 和 时 变 路 网 作 
为 为 衡量 因素 建立 了 卡车 与 无 人 机 联合 配送 的 路 径 优化 模型 ， 
引入 k-means 聚 类 , 设计 混合 粒子 群 算法 对 算 例 进行 仿真 分 析 。 


1 题 描 述 与 模型 建立 


1.1 问题 描述 
本 文 研究 的 问题 


体 可 描述 为 某 物流 企业 有 一 个 配送 中 
心 ， 为 节约 配送 成 本 和 提高 配送 效率 ， 采 用 卡车 与 无 人 机 作 
为 配送 工具 ， 为 所 属 企业 的 客户 进行 配送 。 卡 车 从 配送 中 心 
搭载 足够 数量 的 无 人 机 出 发 ， 完 成 所 属 卡 车 和 无 人 机 客户 的 
配送 任务 后 返回 配送 中 心 。 无 人 机 必须 从 客户 点 或 者 配送 中 
心 完成 起 飞 和 降落 。 根 据 客户 的 位 置信 息 ， 采 用 k-means 聚 
类 形成 卡车 的 初始 配送 路 径 ， 最 后 根据 无 人 机 的 载重 和 飞行 
距离 的 限制 ， 进 行 无 人 机 客户 的 安排 。 无 人 机 在 单 次 飞行 中 
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的 时 间 ; 这: 及 卡车 到 达 客 户 点 i 的 时 间 ; ww: 大 无 人 机 到 
达 客 户 点 i 的 时 间 ; ”57 :卡车 服务 客户 点 集合 ，s :无 人 机 服 
务 客户 点 集合 ;5s : 无 人 机 站 点 的 集合 : 4d;: 客户 i 的 货物 
需求 ， Vr : 卡车 的 行驶 速度 ; 丸 : 单位 产品 的 价值 2 : 单 
位 产品 的 利润 5;: 客户 点 i 的 总 价值 ，Vi;: 客户 点 i 带 来 
的 潜在 价值 ， 8 : a 7 :潜在 客户 
数量 ;so : 信息 扩散 强度 : 信息 扩散 深度 4 : 影响 规 
We ee 


罚 


机 


罚 成 本 系数 及 延迟 惩罚 成 本 系数 ， s,0 : 


分 别 为 卡车 和 无 人 


的 等 待 成 本 系数 。 
3) 决 策 变量 


全 车 辆 从 节点 行驶 到 节点 j 
"10, 否则 
人 无 人 机 7 从 节点 闵行 到 节点 j 
“【0， 否 则 


1.3 客户 价值 及 时 变 路 网 因素 分 析 
1.3.1 客户 价值 分 析 


可 以 完成 多 个 客户 的 配送 ， 超 出 无 人 机 载重 距离 限制 的 客户 
卡车 完成 ， 同 一 路 径 的 客户 由 卡车 和 无 人 机 协同 完成 ， 示 
图 如 图 1 所 示 。 


训 


OO 客户 


配送 中 心 


一 一 > 卡车 行驶 路 径 


二 无 人 机 行驶 路 径 


全。 配送 点 


< 一 一 > 无 人 机 是 大 的 飞行 距离 


图 1 卡车 与 无 人 机 联合 配送 路 径 示意 图 
Fig. 1 Schematic diagram oftruck and UAV joint distribution route 
1.2 ”模型 假设 及 符号 说 明 
基本 假设 : 

a) 配 送 中心 和 客户 点 的 位 置 和 需求 量 已 知 ， 配 送 中心 的 
需求 量 为 0; 


而 


™ 


客 


能 
举 
大 
的 


Ei 


式 


b) 卡 车 每 次 只 携带 一 架 无 人 机 ， 以 及 所 属 卡 车 和 无 人 机 
客户 的 包 于 ; 

co) 无 人 机 的 最 大 续航 里 程 和 最 大 载重 已 知 ; 

dj) 无 人 机 每 次 携带 多 个 包 右 ， 可 以 服务 多 个 客户 ， 每 架 时 


无 人 机 可 以 进行 多 次 交付 ; 

e) 无 人 机 和 卡车 保持 匀速 行驶 , 无 人 机 的 速度 快 于 卡车 ; 

人 无 人 机 的 续航 时 间 不 变 ， 不 考虑 无 人 机 飞行 速度 和 载 
重 对 无 人 机 续航 时 间 的 影响 ; 

不 考虑 客户 点 服务 时 间 和 无 人 机 取 货 、 

hb) 卡车 上 有 足够 的 无 人 机 电池 ; 

i) 卡 车 和 无 人 机 不 能 重新 访问 任 
被 服务 一 次 ; 

j) 卡 车 和 无 人 机 均 为 同一 种 型 号 。 

模型 参数 及 相关 变量 如 下 : 
1) 集 合 
Nm={1,2,.… 贡 表示 所 有 客户 点 的 集合 ;N=Nnut{0} 表示 
所 有 路 网 节点 的 集合 ，0 表示 配送 中 心 ; Nn 表示 所 有 卡车 
集合 ;MN 表示 所 有 无 人 机 集合 。 

2) 模 型 参数 

2 : 卡车 的 最 大 载重 ， Co : 无 人 机 的 最 大 载重 ，z : 
无 人 机 的 最 远 飞行 距离 ， f1: 卡车 单位 固定 成 本 ; 户 :无 人 
机 单位 机 身 使 用 成 本 ; b:: 卡车 单位 运输 成 本 ; by: 无 人 机 
单位 运输 成 本 ;， 方 : 
无 人 机 从 客户 i 飞行 到 客户 j 的 飞行 距离 ，#: 


电池 更 换 时 间 ; 


何 客户 ， 每 个 客户 只 能 


到 达 客 户 点 i 


车 从 客户 ;到 客户 > 的 行驶 距离 ， di : 


样 


客户 价值 包括 当前 价值 和 潜在 价值 ， 当 前 价值 与 客户 的 
需求 量 相关 ， 计 算 方法 为 
点 三 二 g qih2 (1) 


jE Nm 
中 ，4% 为 客户 j 的 需求 量 ; 
户 需 求 量 的 比例 系数 。 
客户 的 潜在 价值 主要 与 经 营 能 力 、 企 业 声 誉 和 科技 创新 
力 等 因素 相关 。 配 送 的 效率 及 配送 方案 的 优 务 对 企业 的 声 


为 单位 产品 的 利润 ， 放 1: 为 


党 产生 影响 ， 此 外 当 客户 的 满意 度 高 时 ， 潜 在 客户 的 概率 增 


， 潜 在 客户 数量 增多 ， 客 户 的 潜在 价值 也 随 之 增 大 ， 客 户 
潜在 价值 计算 公式 为 
Y= PI1MGA2 (2) 


，7 为 客户 j 激发 的 潜在 客户 数 


中 ，P 为 客户 i 的 满意 度 


，W=5owy ， 其 中 四 = Po6 ，pr 是 客户 的 重要 程度 ， 所 以 
(2) 可 以 改写 为 式 (3)。 
VY = PiSGoPrit WY qjh2 (3) 
客户 价值 为 当前 的 价值 与 客户 的 潜在 价值 之 和 ， 计 算 公式 为 
S; = > Wi(g; + v) (4) 
JENa 
中 ，W 为 客户 的 权重 ,客户 的 重要 度 越 高 ， 权 重 越 高 ， 这 
可 以 提高 客户 的 价值 ， 推 动 企业 的 可 持续 发 展 。 


1.3.2 时 变 路 网 分 析 


成 
下 


时 变 路 网 下 ， 卡 车 的 行驶 速度 随 交 通 流 变化 ， 卡 车 运输 
本 会 随 着 拥堵 状况 的 增加 而 大 幅度 提高 ， 卡 车 在 拥堵 状况 
运输 时 间 增 加 导致 无 法 满足 配送 客户 时 间 窗 的 需要 。 交 通 


引 
律 


堵 主 要 包含 常 发 性 拥堵 和 偶发 性 拥堵 ， 常 发 性 拥堵 具有 规 
性 ,可 预测 ,在 路 径 配 送 中 可 优化 配送 时 间 段 和 配送 路 径 。 
交通 流 的 变化 一 般 用 明确 的 速度 分 布 函数 表示 ， 本 文 在 


廿 


和 人 研究 的 基础 之 上 ， 假 设 在 松散 时 间 段 时 ， 速 度 分 布 符合 
才 数 正 态 分 布 hv(D) ~N(4,07)， 


在 交通 流 较 大 时 即 一 般 时 间 段 


[高 峰 时 间 段 , 速度 分 布 符合 正 态 分 布 wD ~N(o?) 。 函 数 表 
达 式 为 
1 nv)? el 让 计 
ee 2o ,VElvninymax|,IEthw 
Jo)=] 
e 2 ,YE[vmin,ymax],rE fw,tws 
V27 
QOL,IE 1 Ow,1E tw 
HH=40%tEtW, Ov=10w,tEtw2 (5) 
3,1E 1W3 Ow3,1E tw3 


:202205.00131v1 


chinaXiv 


其 中 ， wwzstw 分 别 表 示 三 个 交通 流 时 间 段 ，9,929; 表示 三 
个 时 间 段 的 速度 期 望 值 ，owo0w,0% 表示 速度 标准 差 。 
1.4 目标 函数 

本 文 考虑 客户 价值 的 卡车 与 无 人 机 联合 配送 时 变 路 径 
标 函 数 为 : 总 成 本 最 小 。 
1.4.1 成 本 函数 的 构成 
本 文 目标 函数 中 的 总 成 本 的 构成 如 下 : 
1) 运 输 成 本 
车 辆 在 行驶 的 过 程 中 会 产生 运输 成 本 ， 无 人 机 在 飞行 的 
过 程 中 会 产生 飞行 成 本 ， 此 部 分 成 本 的 计算 公式 为 

21 -bY mb th Dg (6) 

2) 固 定 成 本 

企业 在 物流 配送 中 ， 会 产生 车 辆 的 发 车 成 本 和 无 人 机 的 
机 身 使 用 成 本 。 计 算 公式 为 


ZZ f12 Dxom 十 f22 >》 yw (7) 


j=0 h=1 j=0 h=l 


3) 配 送 时 间 窗 的 提前 惩罚 成 本 和 延迟 惩 避 成 本 
在 时 变 路 网 的 约束 条 件 下 ， 为 提高 客户 的 满意 度 和 客户 
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2 基于 改进 k-means 聚 类 的 混合 粒子 群 算法 


2.1 改进 k-means 聚 类 算法 

传统 的 k-means 聚 类 算法 ， 在 进行 聚 类 之 前 须 事先 指定 
类 别 数目 k,， 在 配送 车 辆 未 知 的 情况 下 , 无 法 事先 确定 具体 的 
用 车 数量 。 本 文采 用 改进 的 k-means 聚 类 算法 ， 对 客户 点 进 
行 聚 类 ， 采 用 距离 作为 相似 性 评 佑 ， 用 disi(56,5)) 表示 两 个 数 
据 对 象 5 与 5; 之 间 的 欧 氏 距离 ， 计 算 公 式 为 
dist(6,6) = (6 6n) +t (On 5p) (20) 
其 中 ， 为 数据 对 象 总 个 数 。 
2.2 k-means 聚 类 中 心 的 确定 

传统 的 聚 类 中 心 是 随机 生成 的 ， 本 文采 用 位 置 编 码 的 方 
式 对 聚 类 中 心 进行 选取 [19]。 假设 样本 的 数据 规模 为 mw ， 每 个 
数据 对 象 有 ?个 特征 属性 ， 聚 类 个 数 为 x ， 聚 类 中 心 为 
(j=4,2,…,R) ,共生 成 M 个 粒子 , 每 一 个 粒子 的 位 置 是 由 大 个 
聚 类 中 心 组 成 ， 其 位 置 编码 结构 为 
particle(i): location[] = [Xe, Xe2,***, Xek (2 1 ) 
其 中 ， 避 为 第 /类 数据 的 聚 类 中 心 ， 是 一 个 Pp 维 向 量 。 
2.3 算法 设计 
根据 以 上 的 聚 类 算法 ， 可 以 初步 得 出 车 辆 的 使 用 数量 和 


价值 ， 会 出 现 因 提 前 配送 而 出 现 的 提前 惩罚 成 本 ， 因 交通 影 
响 而 延迟 配送 的 惩罚 成 本 。 计 算 公式 为 


Z3= > Faxmax (max (i,17),0)+ Bxmax (max (th,tr )—ba0)] (8) 
1=0 /1 
4) 等 待 成 本 


于 配送 速度 及 路 况 因素 的 不 同 ， 在 配送 的 过 程 中 会 发 
生 卡 车 的 等 待 成 本 和 无 人 机 的 等 待 成 本 。 计 算 公 式 为 


2Z4=ey $Y max(0,wir)+ 0 max (0,wir) (9) 


i=0 h=l 


其 中 ， wi*=tw-t 表示 卡 车 4 在 i 点 等 竺 无 人 机 的 时 间 ， 
wir= 忆 -tir 表示 7 无 人 机 在 i 等 待 卡车 的 时 间 。 
1.4.2 数学 模型 
minZ=ZI+Z2+Z3+Z4 (10) 
of > dixSr+>dixSo<Or (11) 
>， py Xiih <>， bb YUr, VIESs (12) 
iEmhEeEH iEmrER 
> > t+> >》yr=LVIESn (13) 
iEmhEeEH iEmrER 
dix Sp < QD, ViE Nm (14) 
Xin=1,VIEST 
和 045) 
Dyir= > wir=0,vESr,vrER (16) 
ai<ti<bi (17) 
XinE{0,1}, Vi, JE Nm.hEH (18) 
yirE{0,1}, Vi, jENm, rER (19) 


其 中 ， 式 (9) 为 目标 函数 ， 表 示 成 本 最 小 化 ， 第 一 部 分 为 运输 
成 本 ， 第 二 部 分 为 固定 成 本 ， 第 三 部 分 为 配送 时 间 窗 的 提前 
惩罚 成 本 和 延迟 惩罚 成 本 , 第 四 部 分 为 等 待 成 本 ; 约束 式 (10) 
表示 卡车 和 无 人 机 服务 的 所 有 客户 的 需求 量 不 超过 卡车 的 最 
大 载重 ; 约束 式 (11) 表 示 卡 车 完成 无 人 机 起 飞 的 客户 点 后 , 无 
人 机 才 可 以 进行 配送 ;约束 式 (12) 表 示 每 个 客户 只 能 被 无 人 
机 或 卡车 服务 一 次 ; 约束 式 (13) 表 示 无 人 机 配送 客户 需求 量 
在 无 人 机 载重 约束 范围 内 ， 约束 式 (14) 表 示 无 人 机 无 法 访问 
的 客户 由 卡车 进行 访问 ; 约束 式 (15) 表 示 超 出 无 人 机 距离 约 
束 的 无 人 机 无 法 访问 ; 约束 式 (16) 表 示 客 户 的 时 间 窗 约束 ; 约 
束 式 (17)(18) 为 决策 变量 。 


车 辆 所 服务 的 客户 ， 为 进一步 优化 车 辆 的 配送 路 线 ， 引 进 遗 
传 算法 中 的 交叉 和 变异 概念 ， 设 计 混 合 粒子 群 优化 算法 对 模 


2.3.1 混合 粒子 群 优化 算法 
粒子 群 算法 是 一 种 启发 式 优化 算法 ， 模 拟 鸟 群 随机 寻找 


食物 ， 通 过 经 验 和 交流 ， 调 整 自己 搜寻 的 方向 和 速度 ， 从 而 
找到 最 优 解 (71。 
假设 在 一 个 DD 维 的 目标 搜索 空间 中 ,， 有 7 个 粒子 组 成 一 


个 群落 ,其 中 第 i 个 粒子 的 位 置 表 示 为 一 个 D 维 的 向 量 
-CnXiasXo) ; 第 i 个 粒子 的 历史 最 优 位 置 为 
P=(pupe po); 整 个 粒子 群 迄今 为 止 搜索 到 的 最 好 位 置 记 
为 P=(peprs…pp) ; 第 i 个 粒子 的 运动 速度 也 是 一 个 D 维 的 
向 量 Vi-(VoV…Vp) 。 对 于 粒子 i 第 k+1 次 迭代 时 的 粒子 速 
度 和 位 置 表示 为 


w= wo—(wo—w) (22) 


G1= C2 mm (Cum Cn) (23) 


只 =wx 矶 cixrandOx( 万 一 三 )+czxrandOx(ps 一 世 ) (24) 
2 | = 驹 + va (25) 
其 中 ，w 是 惯性 权重 ，“ 和 是 两 个 随机 数 ， 分 别 是 个 体 学 
习 因 子 和 全 局 学 习 因 子 ，rand0 产生 一 个 (0.0) 之 间 的 随机 数 ， 
i 表示 第 i 个 粒子 ， 代表 连 代 的 次 数 ，4 代表 维度 。 
2.3.2 算法 流程 
粒子 群 优化 算法 由 于 本 身 粒子 更 新 速度 的 快慢 ， 会 影响 
全 局 最 优 解 的 产生 ， 为 使 算法 尽快 跳出 局 部 最 优 ， 将 引入 遗 
传 算法 中 的 交叉 和 变异 概念 , 防止 算法 过 早 的 收敛 ,设计 一 种 
混合 粒子 群 算法 (hybrid particle swarm algorithm,HPSA)。 本文 
设计 的 混合 粒子 群 算 法 实现 的 步骤 具体 如 下 : 
Step1: 初 始 化 粒子 群 (粒子 群 共 有 mm 个 粒子 ), 给 每 个 粒子 
赋予 随机 的 初始 位 置 和 速度 。 
Step2: 计 算 适 应 度 值 。 根据 适应 度 函 数 , 计算 每 个 粒子 的 
适应 值 。 
Step3: 求 个 体 最 佳 适 应 值 。 对 每 一 个 粒子 , 将 其 当前 位 置 
的 适应 值 与 其 历史 最 佳 位 置 (zx*) 对 应 的 适应 值 比较 , 更 新 个 
体 最 优 粒子 。 
Step4: 求 群体 最 佳 适应 值 。 对 每 一 个 粒子 , 将 当前 位 置 适 


局 最 佳 位 置 (gw) 适应 值 比较 , 更 新 


用 OX， 对 个 体 粒 子 和 群体 粒子 进行 


用 一 定 概 率 ， 进 行 粒子 的 变异 操作 。 


Step7: 更 新 粒子 的 位 置 和 速度 。 


断 算 法 是 否 结束 。 如 曙 
step2， 若 满足 结束 条 件 ， 输 出 全 
2.3.3 算法 流程 图 


果 未 满足 结束 条 件 ， 返 


混合 粒子 群 算法 流程 如 


到 


开始 


适应 度 值 计算 | 


更 新 个 体 最 优 粒子 和 群体 最 优 粒 子 


Y 
个 体 最 优 0X 交 叉 


v 
群体 最 优 0X 交 叉 


输出 群体 最 优 粒子 


结束 


局 2 ”混合 粒子 群 算法 流程 区 


Flow chart of hybrid particle swarm optimization 


五 


ChinaXiv 合 作 期 刊 


温 廷 新 ， 等 : 考虑 客户 价值 的 卡车 与 无 人 机 联合 配送 时 变 路 径 优 化 方法 研究 第 39 卷 第 10 期 


局 最 优 粒子 。 
优 交叉 。 


2.3.4 无 人 机 路 径 的 局 部 搜索 
通过 改 k-means 聚 类 混合 粒子 群 算 法 ， 将 所 有 的 客户 分 
配给 卡车 ， 得 到 了 最 优 的 配送 方案 。 完 成 所 有 客户 的 分 配 任 
务 后 ， 对 该 客户 再 进行 局 部 的 搜索 形成 无 人 机 的 配送 路 径 。 
无 人 机 路 径 局 部 搜索 步骤 如 下 : 首先 ， 从 所 有 卡车 的 配 
送 路 径 中 , 从 左 至 右 选择 符合 无 人 机 配送 条 件 (最 大 载重 和 最 
大 续航 里 程 ) 的 客户 点 , 然后 在 剩余 的 卡车 路 径 中 选择 两 个 
离 客户 最 近 的 两 个 点 ， 分 别 为 无 人 机 的 发 射 和 降落 点 ， 如 果 
这 个 点 被 其 他 无 人 机 作为 发 射 点 ， 则 从 其 他 客户 中 选择 。 在 
路 径 中 没有 被 无 人 机 选中 的 客户 点 为 卡车 的 服务 对 象 ， 卡 车 
从 左 到 右 为 客户 进行 服务 ， 最 后 分 别 形 成 了 卡车 和 无 人 机 的 
最 优 配送 路 线 。 


3 ” 算 例 实验 及 结果 分 析 


3.1 实验 环境 及 参数 

本 文 的 算 例 测试 ， 在 一 台 Intel(R)Core(TM)i7- 
1065G7CPU, 1.50 GHz, Windows 10(64) 位 的 电脑 上 进行 ， 在 
MATLABR2018b 上 实现 对 混合 粒子 群 算法 的 编程 。 卡 车 的 实 
时 速度 通过 速度 分 布 函数 来 确定 ， 在 求解 参数 中 ，8" 客户 当 
前 价值 阐 值 为 120，V。 客 户 潜在 价值 闵 值 为 20，s? 信 息 传播 
强度 为 0.1，5 信息 扩散 深度 为 0.125， 影响 规模 为 20， 久 4 
单位 产品 的 价值 为 6 元 ， 如 单位 产品 的 利润 为 3 元 ， 车 辆 的 
最 大 载重 为 100 ks ,无 人 机 的 最 大 载重 为 10t8 ， 本 文 设 定 种 
群 N=100, 最 大 迭代 次 数 为 200,w 惯性 权重 为 1，c 个 体 学 
丸子 为 1.5， 全 局 学 因子 为 2。 
3.2 算 例 验证 

以 某 物 流 运 输 企 业 为 例 ( 数 据 从 网 站 Index of 
/chairedistributique/data/ 获 取 )， 要 为 20 个 客户 配送 货物 ， 
配送 中 心 和 客户 的 相关 信息 见 表 1， 编 号 0 为 物流 配送 中 
心 , 编号 1~20 为 客户 点 , 文中 假设 车 辆 的 平均 行驶 速度 50 
kn/h ， 车 辆 的 总 数 为 4 辆 ,根据 车 辆 的 路 线 安排 ， 在 满足 
无 人 机 的 载重 约束 条 件 下 ,合理 安排 卡车 与 无 人 机 的 路 线 
安排 。 


Tr 
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客户 点 相关 信息 
Related information of customer points 

位 置 时 间 窗 重要 度 权重 影响 范围 
40,50 [7,18] 1.0 1.0 0 
45,68 18,20] 1.0 1.0 
45,70 17,19] 0.5 1.0 2 
42,66 [8,9] 0.5 1.0 5 
42,68 17,18] 0.5 1.0 11 
42,65 [7,8] 1.0 1.0 3 
40,69 16,17] 0.3 1.0 15 
40,66 10,12] 1.5 2.0 

38,68 13,15] 1.5 2.0 9 
38,70 13,15] 1.3 1.5 14 
35,66 [8,9] 0.7 1.0 11 
35,69 [9,11] 1.0 1.0 3 
25,85 15,17] 1.6 2.5 13 
22,75 [7,8] 2.0 3.2 7 
20,85 10,12] 1.3 1.2 11 
20,80 19,20] 3.0 4.0 2 
20,85 [7,9] 0.5 0.5 4 
18.75 [9,13] 1 2.8 3 
1S.75: 15,18] 1.3 1.5 13 
15,80 17,19] 1.0 1.0 

30,50 [7,8] 1.3 1.5 4 


对 该 企业 配送 过 程 中 的 交通 状况 进行 调查 ， 一 天 中 具体 
的 交通 状况 见 表 2。 本 文 假设 卡车 的 单位 固定 成 本 为 100 元 ， 
卡车 的 单位 运输 成 本 为 3/ 元 /wn ,无 人 机 的 固定 成 本 为 15 元， 
无 人 机 的 单位 运输 成 本 1 元 /km， 卡 车 的 等 待 成 本 为 3 元 /4， 
无 人 机 的 等 待 成 本 为 1.5 元 /n ， 卡 车 提前 惩罚 成 本 系数 为 1， 

Tab.2 Velocity distribution function 


拥堵 情况 时 间 段 速度 分 布 
畅通 时 间 段 I Inv(1) ~ N(3.8,0.122) 
18: 00-20: 00 
一 般 时 间 段 sy v(t) ~ N(32,0.52) 
14: 00-17: 00 
拥堵 时 间 段 0 v(D) ~ N(15,0.57) 
17: 00-18: 00 
根据 以 上 的 模型 , 利用 MATLAB 进行 编程 , 将 不 同时 间 


段 的 速度 参数 加 入 到 混合 粒子 群 算法 中 ， 并 根据 无 人 机 局 部 
搜索 得 到 卡车 和 无 人 机 的 配送 路 径 ， 见 表 3, 编 程 优化 后 ， 总 
成 本 为 1257.53 元 。 


Ts 表 3 ”路径 规划 结果 

J Tab.3 Path planning results 

Ti 和 - 

CD 编号 卡车 -无 人 机 配送 路 径 成 本 

er 16 ~ 

a 1 0 > 13 > a > 12 > 0 257.85 

OO i 

ED 2 0 > 20 > 17 » 19 w 15 > 0 261.85 
A Bi 

LO) 3 ” 、 159.90 

a 0 > 10 > 1l = 9 > 2 = 1 > 0 

= I 

CN 4 0 的 a a 117.92 

O| 证， 实 线 表 示 卡 熙 运输 , 庶 线 表示 无 人 机 运输 。 

ed 十 比 字 难 一 

GAN 3;21 对 比 实验 


为 验证 动态 路 网 和 客户 价值 对 卡车 与 无 人 机 的 配送 路 径 
的 合理 性 ， 本 文 将 未 考虑 客户 价值 的 静态 模型 与 本 文 考虑 客 
户 价值 的 动态 模型 进行 对 比 实验 ， 将 得 到 的 结果 进行 比较 分 
析 ， 如 表 4 所 示 。 

表 4 未 考虑 客户 价值 静态 路 网 模型 与 本 文 模型 结果 对 比 


Tab.4 Compares the results of static road network model without 


considering customer value with that of this paper 


参数 未 考虑 客户 价值 静态 路 网 模型 本 文 模型 
总 成 本 /元 1150.27 1257.53 
总 客户 价值 /元 2058.473 2404.989 
客户 满意 度 0.7425 0.9005 
根据 表 4 的 对 比 实验 实验 结果 可 知 ， 首 先 ， 在 总 成 本 方 


本 ， 未 考虑 客户 价值 静态 网 络 模型 下 的 卡车 与 无 人 机 的 配 这 
相 比 于 本 文 的 模型 降低 了 9.32%， 减 少 了 物流 公司 的 支出 ， 
提高 了 经 济 效益 。 其 次 ， 从 总 的 客户 价值 和 满意 度 来 看 ， 本 
文 模型 比 考虑 客户 价值 静态 网 络 模型 下 的 卡车 与 无 人 机 的 配 
送 ,总 客户 价值 增加 了 16.83%, 客户 的 满意 度 提高 了 21.28%。 
最 终 ， 由 以 上 的 结果 分 析 可 知 ， 考 虑 客户 价值 和 时 变 路 网 因 
素 的 模型 ， 可 以 更 符合 实际 情况 ， 提 高 公 争 力 ， 实 现 
企业 的 效益 最 大 化 ， 促 进 企业 的 可 持续 发 展 。 

3.2.2 对 比 实验 二 
针对 车 辆 与 无 人 机 的 联合 配送 路 径 问 题 ， 不 同 的 学 者 提 
出 了 不 同 的 算法 进行 求解 ， 如 : 最 短路 算法 与 禁忌 搜索 算法 
相 结合 的 混合 禁忌 搜索 算法 (hybrid tabu search with short-path 
algorithm，SPTS)6I]、 自 适应 大 邻 域 搜索 算法 (adaptive large 
neighborhood search, ALNS)BI、 贪 禁 随 机 自 适 应 搜索 算法 
(greedy randomized adaptive search procedure,GRASP)05 及 混 
合 遗传 算法 (hybrid genetic algorithm, HGA) 08 等 。 但 以 上 的 
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算法 在 不 同 的 方面 均 存 在 一 定 的 局 限 ， 为 验证 本 文 提出 算法 
的 有 效 性 和 合理 性 ， 将 以 上 四 种 方法 对 考虑 客户 价值 的 卡车 
与 无 人 机 的 联合 配送 时 变 路 径 模型 进行 求解 ， 并 将 实验 的 结 
果 进 行 比较 分 析 ， 如 表 5 所 示 。 

表 5 算法 性 能 比较 


Tab.5 Performance comparison of algorithms 


算法 类 别 ”平均 值 最 优 值 。” 计算 时 间 /s ”收敛 迭代 次 数 
SPTS 1528.27 1347.56 70.891 160 
ALNS 1628.53 1548.39 57.294 110 

GRASP 1490.89 1298.36 65.783 80 
HGA 1589.73 1308.61 80.732 120 
HPSA 1473.56 1257.53 60.4785 140 


根据 表 5 的 实验 结果 可 以 直观 的 看 出 , a)SPTS 在 160 代 
完成 算法 的 收敛 ，ALNS 在 110 代 完 成 算法 的 收敛 ， 尽 管 
ALNS 具有 较 快 的 收敛 速度 ， 但 在 配送 成 本 的 最 优 值 方面 比 
SPTS 的 求解 性 能 差 。b)GRASP 和 HPSA 的 实验 结果 相 比 ， 
GRASP 收敛 速度 较 快 , 但 在 计算 的 时 间 上 比 HPSA 略 长 , 平 
均值 和 最 优 值 都 高 于 HPSA 的 求解 结果 。c)HPSA 和 HGA 的 
实验 结果 相 比 ， 在 平均 值 和 最 优 值 的 结果 方面 都 具有 较 强 的 
求解 精度 ,求解 计算 的 时 间 更 快 ,有 利于 更 快 的 获得 最 优 解 。 
综 上 ， 在 与 其 他 算法 的 结果 进行 比较 时 ， 本 文 的 采用 混合 粒 
子 群 算法 在 求解 本 文 的 模型 时 ， 具 有 较 高 的 计算 准确 度 和 较 
快 的 计算 速度 ， 有 效 的 避免 了 算法 在 求解 问题 时 出 现 局 部 最 
优 的 局 限 性 ， 具 有 较 快 的 收敛 速度 和 全 局 搜索 最 优 的 性 
节约 了 物流 企业 配送 的 成 本 ， 提 高 了 企业 的 经 济 效益 。 


4 ”结束 语 


本 文 针对 客户 价值 下 的 卡车 与 无 人 机 联合 配送 的 时 变 路 
径 问题 进行 了 以 下 的 研究 。 

1) 从 配送 模式 方面 ， 本 文采 取 了 卡车 与 无 人 机 联合 配送 
的 方式 ， 将 客户 的 价值 和 路 网 的 时 变性 作为 约束 条 件 ， 建 立 
了 总 成 本 最 小 的 数学 模型 。 将 静态 未 考虑 客户 价值 的 模型 与 
动态 考虑 客户 价值 的 模型 进行 对 比 ， 本 文 的 模型 可 以 更 好 的 
增加 客户 价值 和 提高 客户 满意 度 , 从 而 扩大 企业 的 经 济 效益 。 

2) 为 提高 算法 的 性 能 ， 设 计 了 混合 粒子 群 算法 对 模型 进 
行 求解 ， 将 不 同 的 算法 实验 结果 与 此 算法 进行 对 比分 析 ， 实 
验 结果 可 以 得 出 混合 粒子 群 算法 在 收敛 速度 、 计 算 时 间 和 最 
优 值 的 求解 方面 有 具有 良好 的 性 能 ， 是 解决 卡车 与 无 人 机 联合 
配送 运输 问题 的 一 种 有 效 的 智能 优化 算法 。 

3) 本 文 在 研究 卡车 与 无 人 机 的 路 径 问题 时 ， 只 考虑 了 无 
人 机 载重 因素 的 影响 ， 未 来 的 研究 可 以 将 配送 的 时 间 和 能 量 
的 消耗 进行 综合 的 考虑 ， 同 时 也 可 以 将 天 气 的 不 确定 性 因素 
进行 考虑 ， 结 合 现实 中 的 具体 问题 建立 更 符合 的 模型 。 
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